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ABSTRACT 

Ion  acoustic  waves  have  been  artificially  excited  and  detected  in 
the  Steady  State  Plasma  Facility  at  the  Naval  Postgraduate  School  by  a 
magnetic  excitation  device  and  an  optical  detection  system.   The  dis- 
persion of  these  waves  has  been  studied  as  a  function  of  the  exciting 
frequency  and  the  confining  static  magnetic  field. 

Dispersion  curves  are  presented  which  compare  the  experimental  data 
to  a  new  dispersion  equation  which  is  derived  in  the  present  thesis. 
After  the  validity  of  this  dispersion  relation  is  demonstrated,  the 
measurements  are  used  as  a  diagnostic  tool,  and  five  operating  para- 
meters of  the  plasma  facility  are  deduced  from  the  measurements. 

Theories  are  presented  to  explain  the  behavior  of  the  dispersion 
curves  and  three  of  the  five  operating  parameters. 
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Section  1 
INTRODUCTION 

For  the  purposes  of  the  present  thesis  a  plasma  may  be  defined  as 
a  quasi-neutral  electrically  conducting  mixture  of  positive  ions  and 
electrons  plus  a  number  of  neutral  background  gas  atoms.   Due  to  the 
particular  plasma  production  mechanism  employed  in  the  Steady  State 
Plasma  Facility  of  the  Naval  Postgraduate  School,  the  electrons  are  at 
a  relatively  high  temperature  (22  electron  volts)  on  the  cylindrical 
beam  axis,  while  the  ions  have  thermal  energies  of  the  order  of  .09  ev 
and  the  neutral  gas  atom  energies  are  approximately  .01  ev. 

Due  to  the  high  electrostatic  restoring  forces  created  by  even  a 
small  internal  charge  imbalance,  a  plasma  has  the  ability  to  maintain 
overall  internal  charge  neutrality.   If  a  charge  imbalance  should 
develop  due  to  the  random  thermal  motion  of  the  electrons  and  ions,  the 
particles  will  rearrange  themselves  to  shield  the  remainder  of  the 
plasma  from  the  resulting  electric  field.   Thus  beyond  a  certain  average 
distance  from  the  imbalance  the  electric  field  so  created  will  not  be 
felt.   This  distance  is  called  the  Debye  length  and  is,  of  course,  a 
function  of  the  plasma  parameters.   In  a  true  plasma  the  Debye  length 
must  be  much  less  than  the  smallest  dimension  of  the  plasma  system. 
The  Debye  length  is  given  by:   1 

A  D  -  (  Ne^   } 


where 


€0      is  the  permittivity  of  free  space 
K       is  the  Boltzmann  constant 
is  the  electron  temperature 


N   is   the  equilibrium  plasma  density 

e    is   the  electron  charge 
Ion  acoustic  waves   are   longitudinal   plasma  density  oscillations 
occurring   at   frequencies   well   below  those  of  the  electron  plasma  oscil- 
lations.     The  plasma   ions   and  electrons  move  exactly   in  phase  only  at 
Lu  =  0;    for  ^    /  0    in-phase  motion   is  only   approximately  realized-      For 
the  frequencies   employed   in  the  present    investigation,    however,   the   ap- 
proximation  is   an  excellent  one.      The   ion  oscillations   are   slightly 
larger   in  amplitude  than  those  of  the   electrons.      The  ratio  of    ion  per- 
turbation velocity  to  electron   perturbation  velocity  can  be   shown  to  be 

1   +    ( )        2    ,   where  JO*  is   the  wavenumber.      In   Section   2    it   will   be 

2TT     L  J 

shown  that  ( — =L_£)^,£lO    ,  so  it  may  be  neglected  compared  to  1. 

Due  to  their  small  mass  and  high  mobility  the  electrons  see  these 
low  frequency  waves  as  quasi-static  perturbations.   Ion  acoustic  waves 
are  similar  to  ordinary  sound  waves  in  a  neutral  gas,  except  that  the 
propagation  mechanism  is  the  Coulomb  interaction  between  charged  particles 
rather  than  the  short  range  collision  forces  between  molecules.   Their 
existence  was  first  predicted  by  Tonks  and  Langmuir  in  1929  based  on  the 
fluid  analysis  of  a  two  component  plasma.   3 

Before  proceeding  to  an  account  of  previous  work,  the  Tonks- 
Langmuir  dispersion  relation  for  ion  acoustic  waves  will  be  derived. 
The  following  assumptions  are  made: 

1)  The  plasma  electrons  have  negligible  momenitum  compared  to 
the  plasma  ions. 

2)  The  plasma  ions  have  zero  temperature  and  are  singly 
charged.   Thus  there  is  no  ion  pressure. 

3)  There  is  always  100%  ionization  within  the  plasma. 
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4)  There   is   no   electric   field  associated  with  the  wave. 

5)  The   Debye    length   is   much   smaller  than  the  wavelength, 
that    is    XD«-^-  • 

6)  The   drift  velocity  of  the    ions    is    zero. 

7)  The  unperturbed   plasma  has   uniform  density. 

Under  these  assumptions   the  momentum  transfer  equation  for  the   ions   can 
be  written  as: 
-1 


M 


•V(P  +  p)    =  [-57-  +  V   •^J[  (N  +   n)vj ^-   (N   +   n)VXz 


where 

M  is   the   ion  mass 

p   is   the  electron   pressure  perturbation 

P   is   the  equilibrium  electron   pressure 

V  is   the    ion  velocity    inducted  by  the  wave 

N  is   the   equilibrium   plasma  density  of   charged   particles 

n   is  the   plasma  density  perturbation  of   charged   particles 

e   is  the  electron  charge 
BQ   is   the    static  magnetic   field    strength 

z   is    a  unit  vector   parallel   to   the   direction  of  BQ. 
Subtracting  the  momentum  transfer  equation  for  the    ions    in  the   un- 
perturbed  system  from  the  above  equation  and    linearizing  to  the    limit  of 
small   perturbations   gives: 

-JJ-^P   =   yf  (WO    -ANVgr"    +JlNVr9  (1) 


where  Sl~  is   the   ion  cyclotron  frequency   = 


M 


Now  for  perturbations   of  the  form  e1(-*z   ~     *'  ,    equation   (1)    can 
be   solved  to   give: 
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NVr  =  -Nr  It-  c    9  x    9  ) 

M      >r        il2_^2 
and 

Nve  =  -=£-  2jl.  ( 1 ) 

e  M     Tr        _£2    _o2 

Therefore 

^  •  (Nv9?  -  Wl$„  ^  «^  +  4,|L,  <_I_> 

Taking  the  divergence  of  equation  (1)  gives: 

zl-^2p  =  -ieS^.  (NV)  -  J2SJ-  (NV6r  -  NVr0)  , 


or 


i^2P  =  -i«*j.  (Nv)  -^-<^-  ♦  -g-rE-x   ,  x  rU    (2) 

M  M       ^r^  r      3  r       ^2    _  ^2 

But   by  the   linearized   equation  of  continuity, 

^ .    (nv)  =  -  yy~  =  ic$n 
Therefore  equation    (2)    becomes: 

M  ^r2  r     *5T^    ^2    _^2 

Assuming  there   is   no   radial   electron  temperature  variation, 


and 


T^-     r  tt  *>        "      ja2 -o>2    ar2      r  tt; 


where 


2   =    KT  (definition) 

M 


and  X  i-s  the  electron  temperature. 
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Rearranging   and  defining 


x  =  r 


t#  -^  -  *H  * 


this   equation  becomes: 

x  n' '    +  xn'    +  x2n  =  0 
Its   solution   is: 


iC&z    -A>t),    ,     . 

n   -  nQe  J0(x> 

where  JQ    is   the   zero  order  Bessel   function. 
Therefore   it  must  be  true  that 

A2(-Jrr- 1)(42  -.rjjr?  =  (2-405- ->2 

where  the  boundary  condition  that   the  perturbations   vanish  at  the 
position  r   =  J)    is   used.      (The  use  of  this   type  of   boundary   condition 
will  be  more  fully  discussed   in   Section   2.)     This   can  be  rewritten  as: 

a  2  -  5  79  I  oL 

*       JhSZL    7T       c2 

Now  for  an  infinite  medium  J(-+cO   and  A2  =  — —  or  j£='=^ 


C2 


(3) 


Thus   for  an   infinite  medium  the  application  of   a  static  magnetic   field 
has   no   effect.      Equation   (3)    is   the   well   known  Tonks-Langmuir  dispersion 
relation  for    ion  acoustic   waves. 

From  an  elementary  acoustic  analogy  the  "hot"  electrons  provide  the 
perturbing  and  restoring  pressure  gradients,  while  the  massive  ions  pro- 
vide the   inertia.      This    dispersion  relation  predicts    a  phase  velocity  of 

JfzL.  =   c   =    (  ^T  y&.      Also,    since  -^-     =    £q     ,   the   group  velocity  equals 
the  phase  velocity  and,   there   is   no   dispersion. 

A  similar  dispersion  relation  will   be  derived   in  greater  detail 
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including  additional  effects  in  Section  2. 

Ion  acoustic  waves  can  be  self-excited  or  externally  excited  both 
in  gaseous  discharges  of  low  percentage  ionization  and  in  highly  ionized 
plasmas.   The  existence  of  self -excited  waves  in  low  density  gaseous 
discharges  has  been  observed  by  Alexeff  and  Neidigh  4   and  by  Little  and 
Jones   5   in  the  form  of  standing  waves  in  the  discharge  tubes  employed. 
The  present  investigation  deals  exclusively  with  propagating  waves.   Two 
earlier  experiments  (one  on  standing  waves  and  one  on  propagating  waves) 
will  be  described. 

Ion  acoustic  waves  have  been  excited  by  electrostatic  methods  (grids 
6  and  probes   7  )  and  by  magnetic  perturbations  of  the  local  electron 
diffusion  coefficient   5  .   The  magnetic  perturbation  method  was  used 
in  this  experiment.   Its  theory  of  operation  appears  in  Section  3. 

Once  the  ion  acoustic  waves  have  been  excited,  they  can  be  detected 
by  grids,  probes,  and  light  sensitive  detectors.   In  this  experiment  both 
probe  detection  and  photomultiplier  detection  were  used,  although  only 
the  photomultiplier  detection  system  was  used  to  make  quantitative  mea- 
surements.  Detector  theory  and  operation  are  also  covered  in  Section  3. 

Wavelength  and  attenuation  measurements  were  made  as  functions  of 
exciting  frequency  for  various  static  magnetic  field  strengths.   The 
validity  of  the  dispersion  relation  is  demonstrated,  the  data  are  anal- 
yzed and  the  results  presented  as  dispersion  and  attenuation  curves. 

The  new  dispersion  and  attenuation  theories  derived  in  Section  2 
contain  5  adjustable  parameters.   Comparison  of  the  experimental  data 
with  theory  allows  computation  of  these  5  parameters,  whose  values  are 
presented  as  functions  of  the  static  magnetic  field  strength.   The  ion 
acoustic  waves  are  thus  employed  as  a  diagnostic  tool.   Physical  inter- 
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pretations  of  the  experimental    results    (dispersion  curves,    attenuation 
curves,    and   parameter  behavior)   are  also   given  in  Section  5. 

Wong,   Motley,    and  D'Angelo  have  carried  out   experiments   on  the  dis- 
persion of   ion  acoustic  waves    in  cesium  and   potassium  plasmas    in  strong 
magnetic   fields.      6 

Fried   and  Gould      8     have   shown  that   the  Vlasov  equation  predicts 

ion  acoustic   wave   propagation  even   in  the   absence  of  collisions    in  the 

plasma.      Under  this   condition,    however,   the   analysis   predicts   that   the 

waves   will    be  strongly  damped  due  to   a  "trapping"  of   ions  moving  with 

velocities   close  to   the  wave   phase  velocity.      This   can  occur   if  there   are 

ions  whose  temperatures   are  close  to  the  electron  temperature    (T  =  IT). 

From  equation   (3)    the  wave  phase  velocity   is   c   =    (_5JE_) 2,    and   if  T  =  *f 

M 

there  will  be  ions  meeting  this  condition. 

This  type  of  collisionless  damping  was  first  predicted  by  Landau.   9 
It  was  observed  by  Wong,  Motley,  and  D'Angelo,  and  their  work  provided 
the  first  experimental  proof  of  its  existence.   The  conditions  of  their 
experiment  were  such  that  ion-neutral  collisions  were  negligible  (the 
percentage  ionization  varied  between  40%  and  90%)  .   Their  plasma  was 
produced  by  surface  ionization  of  cesium  and  potassium  atoms  on  a  hot 
tungsten  plate.   This  method  guaranteed  that  electron  and  ion  tempera- 
tures were  approximately  equal.   The  damping  was  so  strong  that  standing 
waves  were  not  observed. 

The  Q-3  Machine  (Princeton)  was  used  for  the  experiment.   The  fea- 
tures essential  to  the  experiment  appear  in  Fig.  1.   The  plasma  was  pro- 
duced at  one  end  of  the  machine,  allowed  to  drift  longitudinally,  and 
terminated  approximately  90  cm  away  by  a  similar  tungsten  plate.   Radial 
confinement  was  achieved  by  adjusting  the  longitudinal  magnetic  field  to 
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12,000  gauss  for  cesium  and  6000  gauss  for  potassium  (these  values  were 
chosen  to  avoid  possible  ion  cyclotron  resonance). 

In  the  Princeton  experiment  an  electrostatic  excitation  method  was 
used,  in  constrast  to  the  magnetic  excitation  method  employed  in  the  pre- 
sent investigation.   The  perturbing  signal  was  applied  to  the  excitation 
grid  and  modulated  the  percentage  ion  transmission  through  the  grid.   The 
movable  receiving  grid  detected  the  density  perturbations. 

Phase  velocity  calculations  were  made  by  measuring  the  slope  of  the 
curve  of  phase  delay  vs.  receiving  grid  position  (the  movable  grid)  for 
the  propagating  wave.   One  result  of  the  experiment  was  that  the  phase 
velocity  was  found  to  be  independent  of  frequency  for  the  range  from  15 
kc  to  100  kc  (ion  cyclotron  frequency  was  140  kc  for  cesium  and  240  kc 
for  potassium).   This  conforms  to  the  predictions  of  the  dispersion  re- 
lation derived  in  the  present  thesis,  although  the  frequency  range  rela- 
tive to  the  ion  cyclotron  frequency  is  different. 

The  damping  observed  by  Wong,  Motley,  and  D'Angelo  was  found  to  be 
in  excellent  agreement  with  the  theory  of  Landau  damping.   In  Section  2 
the  criteria  necessary  for  the  presence  of  Landau  damping  will  be  dis- 
cussed, and  it  will  be  shown  that  it  can  be  ignored  in  the  present  inves- 
tigation. 

Figs.  2  and  3  show  phase  velocities  for  the  ion  acoustic  waves  in 
the  Q-3  experiment  for  propagation  parallel  and  anti-parallel  to  the  plas- 
ma drift  velocity. 

The  drift  velocity  of  the  plasma  was  estimated  from  the  difference 
between  the  phase  velocities  of  the  waves  for  parallel  and  anti-parallel 
propagation.   In  the  present  investigation  the  plasma  drift  velocity  is 
one  of  the  5  adjustable  parameters.   The  magnitude  of  the  drift  velocity 
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is  machine  calculated  from  the  dispersion  data.   Both  parallel  and  anti- 
parallel  propagation  at  1000  gauss  were  studied  in  the  present  investi- 
gation, and  the  difference  in  phase  velocities  for  the  two  directions  of 
propagation  provided  an  independent  determination  of  the  drift  velocity. 
The  method  is  the  same  as  that  employed  by  Wong,  Motley,  and  D'Angelo. 
This  determination  is  discussed  in  Section  5. 

Little  and  Jones  have  made  dispersion  measurements  on  externally 
excited  ion  acoustic  waves  in  the  positive  column  of  a  low  pressure  mer- 
cury arc.   5   The  percentage  ionization  in  the  discharge  was  approximately 
1%,  the  current  through  the  mercury  vapor  discharge  was  fixed  at  9  amps, 
and  the  axial  magnetic  field  strength  was  45  gauss.   Their  excitation 
mechanism  was  the  same  as  that  employed  in  the  present  investigation. 
Essentially  the  magnetic  excitation  device  employed  created  sufficient 
variations  in  the  local  electron  diffusion  coefficient  to  modulate  the 
plasma  density.   (This  method  of  excitation  will  be  explained  in  detail 
in  Section  3.) 

Little  and  Jones  detected  the  waves  by  observing  variations  in  the 
light  output  of  the  plasma  (due  to  the  plasma  density  variations)  by 
means  of  photomultiplier  detection.   This  method  was  also  used  in  the 
present  investigation  and  also  will  be  explained  in  Section  3.   Fig.  4 
shows  a  schematic  diagram  of  the  detection  apparatus  used  by  Little  and 
Jones. 

Their  experimental  dispersion  curves  show  a  dependence  of  phase 
velocity  on  frequency  and  a  low  frequency  cutoff  below  which  propagation 
cannot  occur.   These  results  are  not  directly  comparable  to  those  of 
Wong,  Motley,  and  D'Angelo,  since  the  frequency  ranges  relative  to  the 
ion  cyclotron  frequency  are  entirely  different.   Fig.  5  shows  a  plot  of 
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£j  versus^  for  45  gauss  external  field.   The  frequency  cutoff  condition 
can  be  seen. 

The  results  reported  by  Little  and  Jones  conform  to  the  theory  de- 
rived  in  Section  2,  although  their  agreement  with  their  theory  is  not  as 
good  as  in  the  present  thesis.   This  is  probably  due  to  the  unclear 
nature  of  the  boundary  conditions  in  the  discharge  tube  of  their  experi- 
ment (as  they  have  pointed  out).   The  boundary  conditions  applicable  to 
the  present  investigation  will  be  discussed  in  detail  in  Section  2,  and 
it  will  be  shown  that  they  can  be  unambiguously  specified. 

The  present  investigation  examines  the  dispersion  of  ion  acoustic 
waves  in  a  high  density  plasma  of  high  degree  of  ionization  (approxi- 
mately 90%)  and  considers  the  combined  effects  of  ion  drift,  damping, 
ion  temperature,  and  a  confining  magnetic  field  together  with  finite  beam 
geometry. 

The  present  thesis  has  three  purposes.   First,  an  improved  disper- 
sion relation  is  derived  for  a  cylindrical  plasma  including  the  afore- 
mentioned effects    Second,  this  dispersion  relation  is  tested  against 
experiment.   Third,  from  the  comparison  of  theory  with  experiment,  the 
above  5  parameters  for  the  Steady  State  Plasma  Facility  at  the  Naval  Post- 
graduate School  are  computed  by  least  squares  computer  search  to  give 
the  best  fit  to  the  data.   The  validity  of  the  deduced  values  of  the 
parameters  is  also  demonstrated. 

The  derived  theoretical  solution  to  the  wave  equation  is  a  complex 
equation.   The  real  part  of  the  solution  predicts  the  dispersion  rela- 
tion, that  is  wave  number  as  a  function  of  frequency.   The  imaginary  part 
predicts  the  damping  as  a  function  of  frequency. 

The  solution  to  the  wave  equation  is  derived  by  assuming  the  ion  and 
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electron  pressures  are  scalars.   This  allows  truncation  of  the  set  of 
moments  of  the  Boltzmann  equation  after  the  first  2  moments,  which  are 
the  equation  of  continuity  and  the  momentum  transfer  equation  respective- 
ly.  These  equations  are  combined  with  the  plasma  equation  of  state.   As- 
sumption of  the  azimuthal  and  longitudinal  forms  of  the  wave  behavior 
allows  the  resulting  partial  differential  equation  to  be  reduced  to  an 
ordinary  differential  equation  involving  only  the  radial  behavior.   This 
is  solved  in  series  form  and  yields  the  solution  to  the  dispersion  prob- 
lem when  the  boundary  condition  is  inserted. 

Experimentally  the  waves  are  excited  by  perturbing  the  local  elect- 
ron diffusion  coefficient  perpendicular  to  the  confining  static  magnetic 
field.   Increasing  the  local  electron  cross-field  diffusion  coefficient 
decreased  the  local  electron  density.   Since  the  ions  and  electrons  move 
in  phase,  this  creates  a  deficiency  in  the  local  plasma  density  of  ions 
and  electrons.   This  density  perturbation  propagates  parallel  to  the 
static  magnetic  field.   In  the  present  investigation  the  frequency  of 
the  applied  perturbations  was  varied  from  187  kc  to  352  kc . 

The  wavelength  and  amplitude  of  the  ion  acoustic  waves  are  measured 
as  functions  of  frequency.   The  results  are  reduced  by  machine  computa- 
tion and  fitted  against  theory.   As  has  been  previously  mentioned  the  5 
adjustable  parameters  in  the  theory  are  ion  drift  velocity,  the  effective 
ion  collision  frequency  for  momentum  loss,  ion  thermal  energy,  plasma 
beam  diameter,  and  the  e-folding  distance  for  the  radial  density  profile. 
The  first  3  are  characteristics  of  the  plasma  facility,  while  the  last  2 
are  functions  of  the  applied  static  magnetic  field.   Independent  deter- 
minations of  4  of  these  5  quantities  have  been  made,  and  these  values  are 
compared  to  the  results  of  the  least  squares  adjusted  values  obtained  in 
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the  present  investigation.   The  results  are  presented  as  functions  of  the 
static  magnetic  field,  and  theories  are  presented  to  explain  their  be- 
havior. 
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Section  2 
BASIC  THEORY  AND  DERIVATION  OF  THE  DISPERSION  RELATION 

If  Landau  damping  is  to  be  a  factor  in  wave  propagation,  the  dis- 
persion analysis  must  begin  with  the  consideration  of  the  particle  dis- 
tribution functions,  since  the  mechanism  involved  is  the  resonant  or 
near-resonant  transfer  of  energy  between  the  wave  and  those  particles 
having  thermal  speeds  approximating  the  wave  phase  velocity.   If  there 
are  more  particles  whose  thermal  speeds  are  slightly  less  the  wave  phase 
velocity  than  there  are  particles  whose  thermal  speeds  are  slightly 
greater  than  the  wave  phase  velocity,  there  will  be  more  particles  taking 
energy  from  the  wave  than  there  are  giving  up  energy  to  it.   Then  there 
will  be  a  net  loss  of  wave  energy,  and  Landau  damping  will  exist.   In 
any  system  having  a  Maxwellian  distribution  there  will  be  some  particles 
having  thermal  speeds  approximating  the  wave  phase  velocity,  but  for  the 
particular  plasma  under  consideration  this  is  a  negligible  fraction  of 
the  total  number  of  particles. 

As  will  be  shown  later,  the  phase  velocity  of  these  waves  is  greater 

than  or  equal  to  (  KjL  )  2.   Spectroscopic  measurements  by  R.  Booth  and 

M 

R.  L.  Kelly  of  this  laboratory  have  established  that  ion  energies  are 
less  than  0.5  ev  when  argon  is  used  in  the  plasma  facility.   10   Probe 
measurements  made  by  D.  M.  Gall  have  established  that  electron  energies 
along  the  cylindrical  beam  axis  are  22  ev.   11   Thus  it  may  be  safely 
assumed  that  the  electron  temperature  greatly  exceeds  the  ion  tempera- 
ture.  (This  assumption  is  also  consistent  with  the  results  of  the  data 
reduction  presented  in  Section  4.)   Thus  there  will  be  an  insignificant 
number  of  positive  ions  capable  of  causing  Landau  damping. 

An  estimate  of  the  Landau  damping  may  be  made  by  examining  a  dis- 
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persion   relation  derived   by  Fried   for  the  case   T>^T. 
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where 


and 


m   is  the  electron  mass 


M   is   the   ion  mass. 


For  the  plasma  studied  the  Debye  length  is  of  the  order  of  0.02  mm, 

A2k>2    s>      -14 
whereas  typical  wavelengths  varied  from  5  cm  to  20  cm.   Thus  = — ~—  \>  10 

4  * 
and  may  be  neglected  compared  to  1.   Making  this  approximation, 
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Thus    it   can  be   seen  that   an  upper   limit  on  the  Landau  damping  term    (the 
imaginary  term   in  brackets)    is    (     m  ) 2  =   2.3  x  10      .      This    is   sufficient- 
ly small   that    it  may  be   ignored  compared  to   the  damping   arising  from   ion- 
neutral   collisions.      This   neglect   of  Landau  damping   is   even  more  fully 
justified   if   ion-ion  collisions    (ion  fluid  viscosity)   are  considered. 
The  effect   would  be   to    introduce  disruptive  collisions  which  would  des- 
troy the  phase   relationship  between  the  resonant    ions   and  the  wave.      This 
critical    phase  relationship  could   also   be  destroyed  by   ion-neutral   colli- 
sions,   but   not   appreciably  by   ion-electron  collisions   due  to   the  negli- 
gible electron-to-ion  momentum  transfer.      It   should   also   be  noted   that 
the  thermal   speeds  of  the  electrons   are  much  higher  than  the  wave  phase 
velocities   encountered   in  the  present    investigation.      Thus   the  electrons 
see  these  waves   as   quasi-static   and  do   not   contribute   any  Landau  damp- 
ing.     The  process  of  Landau  damping  therefore  need   not   be  considered   for 
the  purposes  of  this   thesis.      The  derivation  will   neglect  moments   of  the 
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Boltzmann  equation  higher  than  the  second. 

For  the  conditions  in  this  experiment  under  which  ion  acoustic  waves 
propagate,  ions  and  electrons  move  in  phase  due  to  the  strong  electro- 
static restoring  forces,  so  the  instantaneous  ion  and  electron  densities 
are  equal.   Although  the  Coulomb  interaction  is  the  mechanism  by  which 
these  waves  propagate,  there  is  no  macroscopic  electric  field  associated 

with  the  waves.   This  assumption  is  justified  at  any  frequency  if  the 

Bo2 

magnetic  pressure,  ( ),  is  much  greater  than  the  plasma  pressure 

(essentially  NKT).  [l3J  In  this  experiment  ( 2 )>  4000  joules/m3, 

while  NfcY  =  10  joules/m  .   The  following  derivation  will  therefore  omit 
consideration  of  any  electric  field  associated  with  the  waves. 

In  the  following  derivation  the  quantities  W,  BQ ,  T  ,  T  are  assumed 
to  be  functions  neither  of  position  nor  of  time. 
W  is  the  equilibrium  ion  drift  velocity 

v  is  the  effective  ion  collision  frequency  for  momentum  loss 
due  to  ion-neutral  collisions 
These  assumptions  are  justified  in  the  cases  of  BQ  and  T.   They  are 
probably  approximately  true  in  the  case  of  T  as  well.   Since  W  is  defined 
as  an  equilibrium  drift  velocity,  it  must  be  independent  of  time,  and  for 
simplicity  it  is  assumed  to  be  independent  of  position  as  well. 

The  momentum  transfer  equation  for  ions  is  first  linearized  and  then 
combined  with  the  linearized  equation  of  continuity  and  the  plasma  equa- 
tion of  state.   The  usual  form  is  assumed  for  the  azimuthal  and  longitu- 
dinal forms  of  the  wave  solutions,  and  a  series  solution  of  the  resulting 
ordinary  differential  equation  of  the  ion  density  perturbations  is  ob- 
tained.  Insertion  of  the  radial  boundary  condition  yields  the  disper- 
sion relation  which  is  then  compared  to  the  experimental  data. 
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Assuming  the  electron  and  ion  pressures  are  scalar  quantities,  the 
momentum  transfer  equation  for  the  ions  under  the  influence  of  the  per- 
turbation is: 

^  ^(P  +TT  +  P  +  TO  =  [^r  +  (V  +  W)  -^]  [  (N  +  n)  (V  +  W)J     (4) 

+  (N  +  n)  (V  +  W)  ?    -J7(N  +  n)  (V  +  W)X* 
where 

P   is   the  equilibrium  electron  pressure 
M    is   the  equilibrium   ion   pressure 
p   is   the  electron   pressure   perturbation 
If    is   the    ion  pressure  perturbation 
V  is   the    ion  velocity   induced  by  the  wave 
n   is   the   ion  density  perturbation 
and  the   positive   z  direction   is  taken   in  the  direction  of  BQ. 

The  momentum  transfer  equation  for  the   ions    in  the  unperturbed 
plasma  column   is: 

~j  V(P   +TT)    =  [^—  +   W,»Jsj][  NWJ  +  NW-?-_QNWXz  (5) 

Expanding  equations    (4)    and    (5),    subtracting    (5)   from   (4),    and 
assuming  the   density  perturbation   is   small  compared  to   the   unperturbed 
density   allows  the   resulting  equation  to  be   partially   linearized  to   the 
first   order   in  small   quantities. 

-1  n,  »  n        „    c>V      .   77   ^n      .    -    ^n 


+   V    •  ^7    (N   +   n)  (V  +   W) 
+  W  •  V    (NV  +  n V  +   nW) 


+   •?  (NV  +   nW)    -  n.  NVXz 


(6) 
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Since   solutions   of  the   form  e     '  -col)    are   SOught,    the  follow- 

ing operators  may  be   replaced  by  the    listed  variables: 

-^r >    -ico    (  cj  assumed   real) 

-^- —  — •*■     \-M,   {M>  complex    longitudinal   wavenumber) 
^ ^    iM    (  M  complex   azimuthal   wavenumber) 


Also   since  W  is    assumed   to   be   in  the    z  direction,   W  =  Wz   and 

w-V— »w^-. 

The  radial  form  of  the  density  perturbations  is  left  unspecified  for  the 
present. 

Linearizing  the  expression  V  •  V  (N  +  n) (V  +  W)  to  the  first  order 
in  small  quantities  gives: 


V 
where 


V   (N  +  n)  (V  +  W)  =  W  (V^r-  +  -=^-  nV0  +  i'itiV-)  , 


—        \  A        A 

V  =  Vrr  +  vee  +  Vzz. 


Thus  (6)  becomes: 


-1  N    -  v    -  e>N      A 

"jjVCp   +  W)    =    -i«ONV  -    i^n\J  +   WVr  *j~  z  (7) 

+    i^W(NV  +   nW)    +  ?N(V  +   W)    -_n.NVQr    +  _C2  NVr9 
From   (7)    are    immediately   obtained: 

-^-(p  **>«*>-  aw  j»  ^J2(^_  +  .|£-) 

MNVa    = (7a) 

y  (io)_    UW  -?  )2    +il2 

(I*-  +  ^L)(i^-    iiw  -7)    -      J^-SLi?   +  *> 

MNVr    =  — S *£ E (7b) 

(i^-    L&7  -»>)2    +S2 
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N 


NV     = 
z 


■^=-(p   +  '&)    -    iOwn   +   iiw2n   +  >?WN   +  WVr^ 
(  iA  -   iiw  -  P  ) 


(7c) 


Forming  the  vector  difference  of  equations    (7a)    and    (7b)    and  taking 
the  divergence  of  the  resulting  quantity  gives: 


V  •     (NV6r   -   NVr9)    =  -^p 


(1^3-  ij&w  -  v>)2  +.a2 


(8) 


Neglecting  volume   ionization  and  recombination  the  equation  of 
continuity  for  the   ions    is: 

^-    (N  +  n)(V  +  W)    =   -  -£- (N   +  n)   =    itOn 


(9) 


Taking  the  divergence  of  equation   (7)    gives: 


^i-V2(p   +tt)   =   -i«0JJ7 


NV  +  V   •    nW 


J  +  W 


g>N 


] 


+   i^wK?  •    NV  +  ^7  •    nW     +9  V.    NV  +  V  •    nW  J 
-JZ^7-[   NVQr   -  NVrel 


(10) 


Linearizing  equation    (9)    gives: 

V  •    (N  +   n)  (V  +  W)    =  ^7  •    NV  +   ^  •    nW 

Using  equation  (11)  to  evaluate  V  •  (NV)  +  V  •  (NW) 
and  using  equation  (8)  to  evaluate  y  •  (NVgr  -  NVr9), 

equation  (10)  becomes: 


(11) 
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1    ^y2(p   +  If  )   =   -  cd2n  +  cJ^Wn  -   icdJn 

M 

<<V2  -  yjz-)  (P  +V) 


w 


MN 


_a2 


M 


(i<o-  iiW-p)2  +  si2 


(12) 


[  dkxcZ-  iiw  -?)^ 


r  ^z 


]  (P   +1t 


(i«0  -   iJ&W  -  \?  )      +^2: 


Under  the  assumption  of  the  spatial  invariance  of  Y  an&   T  equation 
(12)  becomes: 


(c2   +  C2) 


¥n        l*n       *2    ]f         (i*>-   it*  -?)2  7        *2     ] 

^2  +  rT7"   r2"JL     (-j.   i^.^2    +^2  J    ~^n  J 


,2        _2v       W      g>N 

+  (c    +  c  }  —  TF 


(io-  i&w  -  ?)2  +  si2 


(13) 


+  ij2n  -eA^Wn  -    i*}?n  =  0 

where  c      s     SX       (definition) 

M 


and       C2  =  -£^L     (definition). 

M 


The  differential   equation    (13)  would  be  directly   reducible  to 

Bessel's   equation  of  order  A   if  were  equal   to    zero    (that   is    if  the 

'  *>r 

unperturbed   column  had   a  uniform  radial  density  distribution).      This, 
however,    is   not  the   case   in  the   present    investigation. 

The  mean  free  path  for  the   ion-neutral    collision  process    is   greater 
than   1  m,   while  the   ion  gyro   radius    is    less   than   10        m.      Up  to   several 
ion  gyro   radii  away   from  the  beam  axis  the   analysis   of   Pf irsch  and 
Biermann  predicts   that   for  a  Maxwellian  velocity  distribution  and   a 
Gaussian  distribution  of  the  guiding  centers   of  the  charged  Darticles 
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there  will  be  a  Gaussian  density  profile.   14  Previous  probe  analysis 

by  Gall  and  Oleson  has  fitted  a  Gaussian  radial  density  profile  to  N, 

at  least  to  a  good  approximation.   15  Thus  N  =  NQe~o  r  ,  so  that 

— —  y    =  -2  £  r,  and  equation  (13)  becomes: 
N   c>  r     ^ w 


(c2  +  c2)fP2n    *  J-  lSL  -  *L-  „ir        (i^-  iAw-P)2        )    ,2n) 

i^^Ci^-   i^f-P)    +  ^-i2^h     "J 


-    (c2    +  C2)(2Wr£2) 


(i«i-  iiw  -i?)2  +U22 


(14) 


+^)2n  -O^Wh  -   U3y>n  =  0 
Defining     s   =  X  r 


and 


\ 2  «  r     -a2         -ilf  ^  _  ^(oj-^w  +  l£)  1 


equation    (14)    becomes: 


,2     g^n         _^n_ 


£s2  }s 


O0-J&W  +   i?)2 


s   _  2s3g2lw 

(<0-d,W  +   ip  )    2 


+  n(s2-/<2) 


=   0  (15) 


Assuming  a   solution  of  the  form 

n  =  V^  Z  afi  <ao  *  0) 

Sp      j=o 

and   inserting  this    into  equation    (15)    leads   to  the   recursion  relation 


;j  p   (j    +  J^L  -  4 


l2j       -TOT 


2g 


21"'    2 J  j:   (j  +|/-l):  p(. 


2  .      "      2g 


and 


/>  =  M 


(where  a^  has  been  chosen  as  , — i ). 
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Thus   the   solution  of  equation    (15)    is: 


where 


j=o        2  l/A|  2J  j:   (j  +//W  ):   P   (b) 


n  =  lue^6^8   "^ 


(16) 


and 


b   = 


g   = 


2g 
2<52/&W 


(«0-iw  +   i?)X2 


Equation    (16)   reduces   to   the  Bessel   function  J^  for  the  case  g   =  0 
(uniform  density   in  the  unperturbed   column).      For  this   case  the  disper- 
sion relation  would  assume  the   simple  form   Xx  =  2.405  where  the   boundary 
condition  has   been   imposed  that  r  =  X    is   the  radial   position  at  which  the 
perturbation  vanishes.      The  validity  of  this  type  of   boundary  condition 
was    confirmed  by  the  probe  measurement   discussed   in  detail    in   Section   3. 
Independent  confirmation  that  this   type  of  boundary  condition   is    appro- 
priate can  be  obtained   from  an  examination  of   the   radial   electron   temper- 
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This  curve  is  repro- 


ature  profile  as  obtained  by  Gall  and  Oleson. 
duced  as  Fig.  18  in  Section  5.   It  is  evident  from  an  examination  of 
Fig.  18  that  the  electron  temperature  drops  sharply  (or  at  least  becomes 
small  compared  to  the  value  at  r  =  0)  at  the  position  r  =  7,8  mm.   An 
examination  of  equation  (14)  under  the  condition  that  both  c  and  C  be- 
come small  reveals  that  the  magnitude  of  the  density  perturbation,  n, 
must  also  approach  zero.   The  use  of  this  type  of  boundary  condition  is 
particularly  convenient,  since  no  assumption  need  be  made  concerning 
the  existence  of  either  a  pressure  node  or  a  velocity  node  at  the  plasma 
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boundary.      This   problem  has   plagued    investigators  working  with   low  ion- 
ization gaseous   discharges.      5 

For  the  case   g  /  0  the  dispersion  relation  f or  M  =  0    is: 


00 


3=0      2Jj:  JS  rc--^-) 


(17) 


where   now 


g   = 


2S2^.W 


(s03  +    i^)X2 


and 


\2   =]  -I^-2  _    iT/j*    .  ^(^-£w  t    i?)~ 

A       ~[(^-£W+    i?)2  JL  c2   *  C2 


Equation    (17),   or  more   accurately  the  recursion  relation  between 
the   coefficients,    is  compared   to  the  experimental   data  by  a   least 
squares   direct   search  computer  program  on  a  CDC   1604. 
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Section  3 
EXPERIMENTAL  ARRANGEMENT 

A  schematic  illustration  of  the  Steady  State  Plasma  Facility  at  the 
Naval  Postgraduate  School  appears  in  Fig.  6.   The  plasma  is  produced  by 
a  hollow  cathode  discharge,  which  has  been  investigated  by  Lidsky  et.  al. 
Fl6  The  details  of  its  construction  appear  in  Fig.  7.   The  facility  is 
capable  of  operation  on  helium,  neon,  and  argon,  although  only  argon  was 
employed  in  the  present  experiment.   The  neutral  argon  gas  is  supplied 
through  a  micrometer  valve  assembly  to  the  hollow  cathode.   An  r-f  volt- 
age is  supplied  between  the  cathode  and  the  first  anode  in  order  to  ig- 
nite the  arc.   Once  the  arc  has  been  struck,  the  discharge  is  self- 
sustaining,  and  the  r-f  voltage  is  turned  off. 

The  plasma  so  produced  drifts  through  the  hole  in  the  first  anode 
and  into  the  experimental  tube.  This  section  is  a  3  meter  pyrex  tube  of 
inside  diameter  10  cm.   It  is  axially  aligned  with  the  confining  magnetic 
channel  to  within  1  mm  throughout  the  length  of  the  experimental  region. 

The  optimum  operating  conditions  for  the  facility  are  an  arc  cur- 
rent of  100  amps  and  a  voltage  of  90  volts  between  the  cathode  and  first 
anode.   This  condition  is  most  conveniently  achieved  by  regulation  of 
the  argon  gas  pressure  through  the  micrometer  valve  assembly. 

Differential  pumping  of  neutral  gas  atoms  by  the  diffusion  pumps 
maintains  a  neutral  gas  pressure  of  approximately  3  x  10~5  mm  Hg  through- 
out the  experimental  region. 

The  strength  of  the  confining  magnetic  field  is  continuously  variable 
from  500  gauss  to  10,000  gauss.   Some  experimental  data  were  taken  at  a 
field  strength  of  600  gauss,  but  the  agreement  with  theory  was  quite  poor. 
It  was  suspected  that  some  drift  occurred  in  the  current  supplying  the 
field  coils,  so  that  the  confining  magnetic  field  was  not  constant  with 
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time    (a  necessary   condition   for  taking   a   complete   set   of  dispersion  data). 
These   data  are  not  reported,    and  the  minimum  field   strength  employed  was 
800  gauss. 

Upon  reaching  the  end  of  the  experimental   region,    the  plasma   ions 
and  electrons   strike  the  floating  anode    (so-called  because   its   potential 
is   allowed   to   float)   and   are   recombined.      They  are   then  pumped   away  as 
neutrals  by   the  vacuum  pumps    in  the  burial   chamber.      A  photograph  of  this 
section  appears    in   Fig.    8.      A  side  view  of   the  entire   facility  appears 
in  Fig.    9,    and  a  view  of  the  cathode-first   anode  region  appears    in  Fig. 
10. 
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FIG.  8 
END  VI  2d  OF  PLASMA  FACILITY  SHOWING  FLOATING  ANODS 
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FIG.    9 
SIOE   VIEW    OF   PLASMA   FACILITY 
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FI3.  10 
Z:j   VIS-'  OF  PLASMA  MACHINE  SHOEING  CATHODE  REGION 
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In  order  to  excite  the  wave,  a  short  r-f  coil  was  tightly  wrapped 
around  the  pyrex  tube  as  shown  in  Fig.  6.   It  was  designed  for  optimum 
coupling  to  the  plasma  at  the  center  of  the  band  of  exciting  frequencies 
employed.   This  coupling  effectiveness  was  determined  by  the  length  of 
the  coil,  which  was  approximately  that  of  the  wavelengths  excited.   Ex- 
perimentally it  was  found  that  the  r-f  field  strength  dropped  appreciably 
in  the  last  \  of  the  coil  at  each  end,  so  that  these  outer  regions  were 
ineffective  in  producing  a  perturbing  field.   The  coil  itself  was  approx- 
imately 15  cm  long,  leaving  approximately  7  cm  for  the  effective  field- 
producing  region.   This  was  in  the  center  of  the  band  of  wavelengths 
excited. 

The  coil  was  the  inductive  arm  of  a  parallel  "almost  resonant"  tank 
circuit,  which  was  tuned  (with  approximately  1500  wuf  excess  capacity) 
to  resonance  at  the  applied  transmitter  frequency.   This  slight  departure 
from  the  resonance  condition  was  necessary  so  that  the  transmitter  would 
see  a  load  impedance  similar  to  that  of  the  radio  transmission  antenna 
which  it  was  designed  to  feed. 

The  transmitter  used  to  power  the  tank  circuit  was  a  U.S.  Navy 
type  TAB-7  manufactured  by  the  Westinghouse  Corp.   It  is  capable  of 
producing  a  peak  r-f  power  level  of  2  kw  over  the  frequency  range  of 
100  kc  to  500  kc.   This  full  power  condition  was  not  used  due  to  the 
inability  of  the  r-f  exciting  coil  to  dissipate  the  heat  generated  at 
this  condition.   Typical  operating  conditions  were  a  transmitter  output 
current  of  3-4  amps.   The  Q  value  of  the  tank  circuit  was  estimated  to 
be  about  10  at  the  frequencies  used.   Thus  the  actual  exciting  current 
in  the  coil  was  approximately  30-40  amps.   Direct  measurement  of  the 
r-f  magnetic  fields  produced  yielded  values  of  from  125  to  195  gauss, 
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depending  on  transmitter  power.      As    previously  mentioned,    this   perturbing 
magnetic   field  produced  variations    in  the   local   electron  diffusion  co- 
efficient,   creating  plasma  density  variations. 

Assuming  the  magnetic   perturbation,   B,    is   of   the  form  B   =  B-,sino)T 
and  that  0>_fl  (which  was   true   in  all   cases),   only  the   influence  of  the 
perturbing  field  on  the  electrons   need  be  considered,    since   the  massive 
ions   are  unable  to   respond  directly  to  the  rapidly  varying  magnetic   field, 
while  the  electrons   can   instantaneously  follow   its   variations.      The  dif- 
fusion rate  of  electrons   perpendicular  to  the  magnetic   field    is: 


°B    = 


f    e(B   +  B0)    -1 

L         M        J 


where     Dg   is  the   local  electron  diffusion  coefficient   with  both  B   and 

B0  applied 
D     is   the   local   electron  diffusion  coefficient   with   zero   applied 

magnetic   field 
e      is   the  electron  charge 
m      is  the   electron  mass 

)^     is  the  electron-neutral   collision   frequency 

eB 

Now  under  the  conditions  of  the  experiment,  ( — )  »  1.   Therefore 

m  N 

Dg   =   D  [       -   m    M    -   1    2   or  DgCC   1 5-  . 

"*  |_  e(B    +  BQ)    J  **  (B    +   BQ)2 

e)  ^N  r\ 

Now  for  the  general   diffusion  problem  DB  -^ — ^—  =  (V  (source) 

(under   conditions   of   steady  state).      Replacing  ^     by  — *—  ,   where 

dr2      a1 

a  is  an  average  effective  diffusion  mean  free  path,  there  results 

— L__  =Q    or  N  =   S£.   .   Thus  NoC  ^—  and  (N  +  n)  oC  (B  +  B0)2  or 
a1  DB  DB 
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.              9                                                                                Bisin  cdt    9 
(N  +  n)  oC    (B1sinof+  BQ)^.      Since  B1«BQ,    (N   +  n)    =  N(l   +  — -— )z 


Bo 

26-Lsin  o^  Bj 

=  N(l    + ).      Therefore   n  =   2N  — —  sincOr,    and  the  plasma   den- 

Bo  Bo 

sity  is  modulated  at  the  applied  frequency,  the  strength  of  the  perturba- 
tion being  proportional  to  that  of  the  r-f  magnetic  field.   These  per- 
turbations were  studied  over  the  frequency  range  of  187  kc  to  352  kc  and 
for  static  magnetic  fields  from  800  to  2200  gauss. 

Once  the  waves  have  been  generated,  they  propagate  parallel  to  the 
static  magnetic  field  in  both  directions.   Only  propagation  anti-parallel 
to  the  static  magnetic  field  was  investigated  in  detail.   (One  set  of 
parallel  propagation  measurements  was  made;  this  allowed  an  independent 
determination  of  the  ion  drift  velocity,  W. ) 

Previous  spectroscopic  measurements  have  established  that  virtually 
all  the  light  emitted  by  the  plasma  is  from  excited  neutral  atoms.   10 
Therefore  the  emitted  light  intensity  at  any  point  along  the  plasma  col- 
umn is  directly  dependent  on  the  local  electron  density.   (An  ion- 
electron  light  production  mechanism  would  produce  a  light  intensity  which 
was  proportional  to  the  square  of  the  local  charged  particle  density.) 
Thus  a  photomultiplier  tube  can  be  used  to  detect  the  presence  of  the 
waves.   The  detection  circuit  schematic  is  given  in  Fig.  11.   The  photo- 
multiplier  tube  was  housed  in  a  light-tight  electrostatic  and  magnetic 
shield  manufactured  by  the  Quanta  Corp.   A  1P21  photomultiplier  tube  was 
mounted  inside  this  shield  and  dynode  voltages  were  supplied  by  batteries. 
The  entire  assembly  was  located  in  a  region  of  minimum  magnetic  field 
(less  than  25  gauss).   The  optical  signal  was  fed  to  the  phototube 
through  a  10  ft  flexible  unaligned  optical  fiber  bundle  manufactured  by 
Optics  Technology,  Inc.   This  fiber  bundle  caused  a  61%   loss  in  light 
intensity,  but  the  signal  strength  was  large  enough  that  this  caused  no 
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problem. 

A  Keithley  Model  110  Wideband  Amplifier  with  an  overall  gain  of 
unity  was  used  as  an  impedance  matching  device  between  the  phototube 
and  the  Hewlett-Packard  Model  450A  Wideband  Amplifier  (which  was  set  to 
a  20  db  gain).   In  order  to  improve  discrimination  against  noise  the 
signal  was  then  fed  into  a  Spencer- Kennedy  Electronic  Filter  Model  302, 
which  served  as  a  band  pass  filter  centered  at  the  transmitter  frequency. 

The  next  element  in  the  detector  circuit  was  a  Hewlett-Packard  Model 
400H  Vacuum  Tube  Voltmeter  which  was  used  to  make  amplitude  (attenuation) 
measurements.   Finally  a  Tektronix  Dual  Beam  Oscilloscope  Type  555  was 
used  to  display  the  detected  signal  and  compare  it  with  a  monitor  signal 
taken  directly  from  the  TAB-7  transmitter.   Both  signals  were  subject  to 
various  phase  delays  (due  to  coupling  into  and  out  of  the  plasma,  ele- 
ments in  the  detector  circuit,  etc.),  but  comparison  of  the  phases  of 
the  two  signals  appearing  on  the  oscilloscope  gave  the  relative  phase 
delay  of  the  arriving  signal.   Moving  the  optical  fiber  bundle  along  the 
plasma  column  until  the  phase  of  the  arriving  signal  was  the  same  as  be- 
fore allowed  a  direct  determination  of  the  wavelength  of  the  propagating 
signal. 

Signal  strength  after  amplification  and  filtering  was  of  the  order 
of  100  mv  (depending  on  dynode  voltages  and  transmitter  power).   Signal- 
to-noise  ratios  varied  widely  depending  primarily  on  the  electronic 
filter  settings  and  on  the  transmitter  power.   Individual  filter  settings 
at  each  frequency  employed  were  necessary  in  order  to  discriminate 
against  the  noise  background  in  the  plasma.   The  frequencies  selected 
for  study  were  chosen  in  the  quieter  regions  of  the  plasma  noise  spec- 
trum.  There  was  a  marked  variation  in  the  location  of  these  regions 
depending  on  the  static  magnetic  field  strength.   A  trial  and  error 
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approach  proved  to  be  the  most   productive  method   for   finding  their   loca- 
tion. 

A  typical  oscilloscope  photograph   is   shown   in  Fig.    12.      The   upper 
trace  shows   the  detected  signal    (after  amplification  and   filtering). 
Part   of   its   amplitude   is  due   to   an  unavoidable   inductive   pickup  condi- 
tion.     This   part   was,    of  course,    subtracted   from  the  raw  data  before 
analysis.      Some  plasma  noise   can   be  clearly   seen  within  the  signal    en- 
velope.     The   signal-to-noise  ratio   here    is   slightly  greater  than   1.      The 
lower  trace   is   a  comparison  signal   taken  from  the  TAB-7   transmitter. 
This   photograph  was   taken  for   a  250   kc   signal.      Deflection  sensitivity 
of  both  traces   was    100  ow/cm,    and  sweep   speed  was   5  iAsec/cm. 

Wavelength   and   amplitude  measurements   were  taken  as    functions  of 
frequency  for  various  magnetic   field  strengths.      The  results   provide 
the  raw  data  from  which  the  dispersion   and   attenuation  curves   presented 
in  the   next   section  were  derived. 

To    insure  that  modes   other  than  theM=  0  were   not  excited,    probe 
measurements   were  made   in  an   attempt  to   detect   the   presence  of  the 
M  =   1,    2,    3  modes.       (The  form  of  these  modes    is   sketched  below.) 


h~-\ 


yu=2 


/A.-3 


At  maximum  detection  sensitivity  no  such  waves  were  detected.   Since  the 
excitation  mechanism  had  azimuthal  symmetry,  modes  other  than  u=   0 
should  not  be  present.   Actual  confirmation  that  the  detected  wave  was 
actually  the  U  =  0  mode  was  made  by  probe  measurement.   The  perturbation 


48 


Figure  12,  Typical  Oscilloscooe  Trac 


ce 
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was  observed  to  have  a  maximum  at  r  =  0  (on  the  beam  axis)  and  was  ob- 
served to  decrease  as  the  probe  was  moved  out  from  the  axis.   This  be- 
havior would  only  be  possible  for  the  /*=  0   mode.   All  other  modes  would 
exhibit  a  null  at  r  =  0.   This  probe  measurement  also  forms  the  basis  for 
the  comparison  theory  for  the  cutoff  distance  (that  is  Ji) .      These  mea- 
surements showed  no  detectable  perturbations  beyond  a  distance  of  15  mm 
from  the  beam  axis. 

No  standing  waves  were  observed  to  exist.   This  was  probably  due 
to  the  effect  of  the  burial  chamber  surrounding  the  floating  anode.   In 
the  burial  chamber  ions  and  electrons  are  recombined  and  the  resulting 
neutral  gas  pumped  away.   Thus  the  plasma  charge  density  is  low,  and 
the  neutral  gas  density  is  high.   Both  of  these  conditions  would  tend 
to  preclude  the  establishment  of  a  standing  wave  pattern.   Reflection 
of  the  waves  from  the  floating  anode  could  produce  a  standing  wave  pat- 


tern, but  this  effect  has  never  been  observed. 
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Section  4 
EXPERIMENTAL  RESULTS 

Given  a  fixed  static  magnetic  field  strength,  a  set  of  data  consists 
of  <X  and  G   as  functions  of  cu.      The  raw  data  are  reduced  in  the  follow- 
ing way.   Assuming  cO  is  real,  let  Mi    -  &(.  +   i.  A    (  oc  and  jQ    both  real). 
First  the  measured  wavelengths  are  reduced  to  dimensionless  variables 
by  defining  y  =      .   Knowing  the  wavelength  allows  direct  calculation 
of  the  damping  decrement,  A    .   If  two  intensity  measurements  are  made 

one  wavelength  apart,  =  q,     r\    1  ol  where  |z,  -  z  I  =  '"  . .  .   Thus 

In  & 

I,      _  W* 


—  =  e    *  ,  from  which  A   can  be  calculated  directly.   The  damping 


lo 

decrement  is  then  reduced  to  a  dimensionless  variable  by  defining 

z  =  iV.   .   The  frequency  is  reduced  to  a  dimensionless  variable  by  de- 
fining x  =  -j£-   .   For  a  given  magnetic  field  strength  the  entire  set  of 
data  (that  is  y  and  z  vs.  x)  are  read  into  a  direct  search  computation 
program.   The  5  adjustable  parameters  W,  v?  ,  C,J<,  and  S    are  reduced  to 
dimensionless  variables  by  defining 

<->»    W      -/  _   y?     „_  _  _C_     s  .       c      ,     Sc 

The  computation  program  stmultaeeously  adjusts  these  parameters  to 
find  the  dispersion  function  yielding  the  least  square  error  when  com- 
pared to  the  data.  In  addition  to  performing  this  calculation,  graph- 
ical displays  are  also  presented  for  y  vs.  x  and  z  vs .  x.  The  plotted 
squares  represent  the  data  points,  and  the  continuous  curve  represents 
the  dispersion  relation  using  the  5  parameters  just  computed.  Results 
are  presented  from  800  to  2200  gauss  in  the  following  curves,  the  para- 
meter being  BQ(web/m2). 
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The  following  4  curves  present  the  bahavior  of  If    (the  dimension- 
less  drift  velocity)  and  Q"  (the  dimensionless  ion  thermal  speed)  plus 
the  non-normalized  parameters  W  and  C  as  functions  of  the  magnetic  field 
strength.   Independent  measurement  of  W  was  made  at  1000  gauss  by  study- 
ing parallel  propagation,  and  the  results  of  the  data  reduction  for  that 
set  of  dispersion  data  gave  a  value  within  12%  of  that  obtained  from  the 
anti-parallel  analysis.   This  will  be  discussed  more  thoroughly  in  Sec- 
tion 5.   However,  no  theories  are  presented  to  explain  the  behavior  of 
either  parameter.   A  program  was  prepared  to  fit  a  least  squares  straight 
line  to  the  deduced  parameters  as  a  first  approximation.   These  lines 
are  the  continuous  curves,  and  the  squares  represent  the  fitted  para- 
meters for  each  magnetic  field  strength.   The  standard  deviations  of 
both  slope  and  intercept  for  these  plots  appear  in  Table  1.   The  units 
of  each  of  the  non-normalized  quantities  appear  on  the  individual  curves. 


Table  1 


Quant  ity 
std.  dev.  of  slope 
std.  dev.  of  intercept 
slope 
intercept 


r^ 

c 

W(m/s) 

ion 
energy  (ev) 

.00168 

.0741 

12.2 

.120 

.000264 

.0116 

1.9 

.0188 

.00109 

.415 

7.92 

.727 

-.0124 

.0614 

-90.0 

.00635 
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The  next    set   of   6  curves   represents    the  behavior  of  3    (the  dimen- 
sionless   effective    ion  collision   frequency   for  momentum   loss),     §,    (the 
dimensionless   radial   cutoff  distance  for  the  perturbations),    and  €     (the 
dimensionless  e-folding  parameter  for  the   unperturbed   radial   plasma  den- 
sity distribution)    plus   the  corresponding  non-normalized   parameters  */ , 
J^ ,    and    §  .      Again  the   squares   represent   the   data  points   deduced   from  the 
data   reduction,    and   the  continuous   curves   represent   the   theories    for  the 
behavior  of  these  parameters.      These  theories   are  presented    in  detail    in 
Section  5. 

...  _  S  ' 

By   definition        "j       e  — *-— 

and  6     =  -^-       . 

_f7_ 
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The  following  set  of  curves  depects  the  behavior  of  the  dimension- 
less  phase  velocity,  _ £j_  ,  and  the  dimensionless  group  velocity, 

xtt—  >  for  each  value  of  static  magnetic  field  strength.   Again  the 

c  &A' 

parameter  is  the  value  of  B   (web/m^).   The  squares  represent  the  normal- 
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ized  data  points  and  the  continuous  curves  represent  the  values  of 
and  •?£      as  taken  from  the  dispersion  relation  using  the  values  of  the 
5  parameters  appropriate  to  the  particular  value  of  BQ. 

The  modulation  feature  of  the  TAB-7  transmitter  was  inoperative, 
so  no  group  velocity  measurements  were  made,  and  there  is  no  data  with 
which  to  test  the  theory.   For  this  reason  no  data  points  appear  on  the 
group  velocity  curves. 
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The  last  4  curves  represent  the  general  behavior  of  the  dispersion 
function  throughout  the  first  quadrant  of  X  ,  <Ospace.   Plots  are  pre- 
sented for  y  =  — —  vs.  x,  z  =  (  c   vs.  x, u-  vs.  x,  and  _L_  £„ 

vs.  x.   The  1000  gauss  magnetic  field  setting  was  typical,  and  these 
plots  represent  the  dispersion  function  with  the  5  parameters  having 
values  appropriate  to  that  setting.   Since  the  1000  gauss  data  points 
have  already  been  presented  in  the  previous  curves,  the  data  are  omitted 
from  these  4  curves. 

The  first  3  curves  are  for  the  complete  quadrant,  but  the  group 
velocity  is  depicted  only  for  the  regions  in  which  the  damping  is  small, 
that  is  where  ion  acoustic  waves  might  be  useful  for  communication  pur- 
poses.  These  are  also  the  regions  where  rg  >  0.   In  the  regions  of 
anomalous  dispersion  (.-^n      4*  0)  tne  concept  of  group  velocity  becomes 
vague,  since  it  no  longer  necessarily  represents  the  velocity  of  propa- 
gation of  a  signal.   18 
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Section   5 
INTERPRETATION  OF  EXPERIMENTAL  RESULTS 

Before  proceeding  to   a  more   detailed   analysis   of   the   experimental 
results,   the  basic   behavior   of   the  dispersion  curves    and  the  attenuation 
curves   will   be  discussed.      Consider   first   Fig.    13;    here   a  non-dimension- 
less   plot  of    0\    vs.   6«3  is   presented.      For   simplicity,    and   in  order  to 
isolate  the  effects   of  the   various   parameters  on  the  basic  dispersion 
curve,   consider   first   the  case  of  W  =   0    (zero   ion  drift   velocity); 
P=   0    (no    ion-neutral   collisions   and^jQ  =   0<   +   i#  is    real);    C   =   0    (cold 
ions)  ;    and  5   =   0    (uniform  radial   density  distribution    iij  the   unperturbed 
plasma   column);    the  dispersion  relation,   equation    (17)    reduces   to: 


_  .  i)  <«    -  -7-)   -  r-jj,  , 


where  o(   is  the  first  root  of  JL,  N   =  2.405... 

O  O  '   o 

This  equation  is  plotted  in  Fig.  13  for  the  4  distinct  cases: 
XL  =  0 ,  X  -  CO         (zero  magnetic  field,  infinite  medium) 
SX.  ^   0,  Jl  =  CO  (non-zero  magnetic  field,  infinite  medium) 

J*L   =0,  x  ^  co  (zero  magnetic  field,  finite  medium) 

_/L  t*  0,  £  ^   aO    (non-zero  magnetic  field,  finite  medium) 
The  cases  presented  in  Fig.  13  are  elementary  and  are  presented 
as  an  introductory  clarification  for  the  more  detailed  analysis  which 
follows . 

From  Fig.  13  it  is  evident  that  the  presence  or  absence  of  a  mag- 
netic field  has  no  effect  on  the  dispersion  in  an  infinite  medium.   The 
reason  is  that  the  perturbational  velocities  Vq  and  V  are  always  zero 
for  this  case  (see  equations  4a  and  4b).   Only  Vz  is  non-zero,  and  since 
it  is  not  a  function  of  SL   (see  equation  4c),  there  is  no  dispersion. 
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The  term  cutoff  will   be   used  to  describe  the   propagation  condition 
under  which  the   phase  velocity  of  the  wave   is    infinite    (or  when  the  real 
part   of   the   propagation  constant,   o(  ,    is    zero).      The  term  resonance  will 
be   used  to  describe  the  condition  under  which  the  phase  velocity   is    zero 
(  0(  is    infinite)  . 

For   a  finite  medium  the  high  frequencies    (and   hence  small  wave- 
lengths)   can  easily  propagate   down  the   plasma  column.      As  the  wavelength 
increases   there    is   a  point    at  which  the   plasma  can  no    longer  transmit 
the   slgnal.      «,   appears   as    a  eUto«   at   f„   ^   .    «£_   ,    ana 
means   that   at   frequencies    less   than  O , ,  o(     <0  meaning   propagation   is 
impossible.      This    is   a  consequence  of  the   location  of   the  boundary   sur- 
face   (the  cylinder  of   radius^). 

Finally   for  the  case  of   a   finite  medium  with   a  non-zero  magnetic 

field,    there   is    a  resonance   at   the   ion  cyclotron  frequency.      At   this 

frequency  all   the   ion  momemtum  appears   as   cyclotron  motion.      Again 

2    2 
there    is   a  cutoff,    but  this   time  at  u)  n   =    (Si?   +  )    •      At    fre- 

2  A2 

quencies«CL<(€^C^2  the  plasma  beam  cannot  transmit  the  signals,  and  prop- 
agation  is  impossible  since  o(  s  P« 

In  the  following  discussion  it  will  be  assumed  that  the  dispersion 
relation  is  exactly  the  Bessel  function  of  order  zero.   This  is  equiva- 
lent to  the  assumption  that  the  unperturbed  plasma  column  has  a  uniform 
radial  density  profile.   This  will  considerably  simplify  the  discussion, 
since  the  infinite  series  for  the  dispersion  relation  need  not  be  consid- 
ered.  For  this  case  the  relation  \X   =  2.405  holds.   In  actual  practice 
the  perturbations  differ  from  pure  Bessel  functions  by  no  more  than  about 
10%,  so  the  approximation  is  justified  for  discussion  purposes. 

From  this  point  on  the  assumption  will  be  made  that  the  plasma  is 
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finite  in  the  radial  direction  (^  ?   oo  )  and  that  there  is  a  non-zero 
magnetic  field  (-X1  /  0).   Also  all  variables  will  be  normalized  (in- 
cluding the  5  parameters)  to  the  dimensionless  values.   These  definitions 
are  repeated  here  for  clarification. 


X     =       * 

'V-a- 

-   -    /2>c 

?  = 

w 

-CL 

-C2. 

c 

1  =  ? 

c 

^"-a^ 

e  = 

cS 
J2. 

The  effect  of  ■#  (ion  drift)  is  easily  understood,  since  it  intro- 
duces a  Doppler  shift  due  to  the  moving  medium.   In  this  experiment  fi 
and  Oj  were  anti-parallel,  so  the  phase  velocities  are  slightly  reduced 
by  the  drift.   This  is  illustrated  in  Fig.  14. 

A  set  of  measurements  was  made  with  yi    and  <-i  parallel  at  a  field 
strength  of  1000  gauss.   Direct  measurement  of  the  difference  in  phase 
velocities  for  the  parallel  and  anti-parallel  cases  (not  the  use  of  com- 
puter data  reduction)  gave  an  entirely  independent  determination  of  the 
value  of  the  ion  drift  velocity.   This  value  was  within  12%  of  that  ob- 
tained by  the  computer  data  reduction  process,  and  this  agreement  lends 
strong  confirmation  to  the  values  indirectly  deduced  from  the  measure- 
ments . 

The  effect  of  ion  temperature  is  also  easily  understood.   Since  ion 

acoustic  waves  are  density  perturbations,  they  propagate  at  a  velocity 

proportional  to  (  pressure  term  ^   passing  from  the  case  of  cold  ions 

inertia  term 

to  warm  ions  will  increase  the  pressure  term  by  a  small  amount  equal  to 
the  ion  pressure.   Thus  the  phase  velocity  will  be  higher  for  a  given 
frequency;  that  is  M  will  be  reduced.   This  effect  is  illustrated  in 
Fig.  15.   The  solid  curve  represents  the  case  of  cold  ions,  and  the 
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dashed  curve  represents  the  case  of  warm  ions. 

The  effect  of  T  is  more  involved.   If  ion-neutral  collisions  are 
considered,  there  will  be  a  decrease  in  the  wave  phase  velocity  to  a  value 
smaller  than  in  the  collisionless  case.  J  171  In  the  present  investigation, 
however,  this  effect  is  of  the  second  order.   The  more  important  effect 
of  collisions,  however,  is  that  a  damping  mechanism  has  been  introduced. 
This  means  that  there  is  no  longer  an  undamped  resonance  at  x  *  1.   Rather 
the  resonance  is  diffuse,  and  the  maximum  wavenumber  is  limited  due  to 
the  damping.   Waves  can  also  exist  in  the  previously  forbidden  region, 
1  <  x  <  — — —  ,  although  here  they  are  heavily  damped.   This  effect  is 

illustrated  in  Fig.  16.   Although  attempts  to  excite  and  detect  waves  in 

c*5  9 
the  region  1<  x     z  were  made,  the  heavy  damping  plus  the  long  wave- 
lengths in  this  region  precluded  their  detection. 
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Ion  acoustic  waves    are  produced  by  the   Coulomb    interaction    (a   long 
range    interaction).      Propagation   is   possible  only  when  the    ion-neutral 
collision  frequency   is   small   compared  to   the  wave  frequency.      The   few 
collisions   present  disrupt   the   long-range  ordering,    and  damping   increases 
with   increasing   collision  frequency.      Typical    ion-neutral   collision  fre- 
quencies   for  this    investigation  were   less   than  500   sec"    ,    so  the  conti- 
tion  of  heavy  damping  would   be  encountered  only  at   frequencies   much 
less   than  those   studied    (plus   the  heavily  damped  waves    in  the  region 

Due  to   the   apparent    invar iance  of   "jtf    with  magnatic   field,    it    is   sus- 
pected that  the  mechanism  producing  the   drift    is   an  electric   field    in  the 
interior  of  the  plasma,    rather   than   a  fluid  pressure  differential   between 
anodes.      If  the   ion  drift   were  due   to   the    latter  mechanism,    its  value 
would   depend  on  the  magnetic   field  strength,    since  the   neutral   gas    in- 
jection pressure  varied  with  the  field   strength.      Unfortunately  previous 
electric   field  measurements   were  not   available,    so    a  confirmation  of  the 
mechanism  which  produces   the   drift   cannot   be  made. 

For   ion-neutral   collisions,    A.        =  —     — t ^  I 

lon     fTNg^-c        o 

where  A .     is  the  ion  mean  free  path 
ion 

N_.     is  the  neutral  gas  density 
CT"*     is  the  collision  cross  section  for  the  process 

(V  )  2     is  the  mean  ion  thermal  speed 

Knowing  0~"  and  N   (which  is  manually  adjusted  for  each  magnetic 
c     g 

field  strength  to  give  optimum  machine  operation),  a  straightforward 
calculation  gives  0     =  CT*P  x  246  x  10   where  P   is  the  neutral  gas 

o  o 

pressure   and  ^7"    =  . 
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Thus     =  41.7  x  10' 


Bo 


In  the  above  calculation  it  is  assumed  that  there  is  a  100%  excess  ion 
momentum  loss  per  collision.  This  assumption  is  not  strictly  true,  but 
it  is  certainly  reasonable,  since  a  large  fraction  of  the  directed  mo- 
mentum is  lost  in  any  collision.  It  is  also  assumed  that  the  ion  per- 
turbation velocities  are  small  compared  to  the  mean  ion  thermal  speed. 
Thus  the  collision  frequencies  will  not  depend  on  the  strength  of  the 
perturbation. 

The  above  equation  appears  as  the  theory  on  the  ion  temperature 
curves  in  Section  4.   The  value  of  CT"~  used  was  taken  from  the  data  re- 
duction program  appropriate  to  the  magnetic  field  strength. 

Since  ST\K   1  (the  ions  are  much  colder  than  the  electrons),  the  dis- 
persion equation  is  quite  insensitive  to  changes  in^~  .   This  can  be 
seen  by  examining  the  data  points  of  the  CT"  curves  in  Section  4.   It  is 
considered  by  the  author  that  the  deduced  values  of  ^T"  (and  hence  the 
ion  energies)  should  not  be  interpreted  too  literally.   In  the  C"*  vs.  B 
curve  the  deduced  data  points  are  rather  widely  scattered.   This  inac- 
curacy is  reflected  not  only   in  the  ^T*  vs.  B  curves,  but  also  in  the 
"^  vs.  BQ  curves,  since  4°<^T.   But  it  can  be  seen,  at  least  qualita- 
tively, the  observed  damping  results  primarily  from  ion-neutral  colli- 
sions. 

Variations  in  B   affect  the  propagation  characteristics  (that  is 

2  2 
(X     i 

C&2   =    (H^  +  j-k — )^>  since  jfL  is  of  the  same  order  of  magnitude  as 

2 .   An  analysis  by  Wong  119   concludes  that  the  transverse  boundary 

conditions  are  important  only  if  the  ion  cyclotron  radius,  P  .,  is  com- 
parable to  the  plasma  beam  radius.   At  1000  gauss  P.    =  3.2  mm  and  J(    -  15 
mm,  so  this  criterion  is  approximately  met. 
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Neglecting     &    (uniform  unperturbed   column)   Wong  has   shown  that 
1         _4*2_   ?'f^    ,    ?  T  47t2      ,     <Xo2        ^AV:n 

1  "  FV  Z  ^      2  *  '  ?>?     i^F     iW"  ' 

where  Z      is   the  derivative  of  the   plasma  dispersion  function 
Z(jp    S   exp(-^2)^   ifiP     -    2    f  exp^^j 

with  respect   to    its   argument. 

In  the   present    investigation,  4i    *°»     </       °  CS — //      2T      . 

^7T2        S     4*2/2       ^    ^T 

Thus   Z*    (  ■  ■  ■  )    =  2-L-  +    ^°  L L — ,  and  the  effect  of  the   radial   boundary 
>&V  X  Jl   2       '  ' 

is    important  only   if   2— i—  =  — °  /  L —  .      Since  this    is    indeed  the  case, 

f  ^  2  " 

.   (X  Q2a  .2  m 
it   must  be  true   that   T  =  — ^LJ: — !_  ,   or  T  =   2  ev. 

2i2 

This  is  a  lower  limit  for  the  ion  temperature,  since  any  tempera- 
ture appreciably  lower  would  predict  that  43  ^   is  not  a  function  of  BQ. 
This  estimate  is  considerably  higher  than  the  results  of  the  data  re- 
duction (.05  ev  to  .17  ev),  and  the  spectroscopic  measurements  of  Booth 
and  Kelly.   10  For  that  reason  it  must  be  stated  that  the  ion  energy 
cannot  be  determined  accurately  by  this  experimental  method.   It  is 
probable,  however,  that  the  ion  energies  lie  somewhere  between  .05  ev 
and  2  ev.   The  ion  temperature  has  been  found  to  increase  with  distance 
from  the  plasma  source  in  a  helium  plasma,  varying  in  value  from  about 
1  ev  at  the  first  anode  to  about  8  ev  several  meters  downstream.   20 

The  increase  of  ion  energy  with  increasing  magnetic  field  is  be- 
lieved to  be  due  to  the  increased  electron- ion  collision  frequency  in 
the  denser  beams  resulting  from  higher  field  strengths. 

It  should  be  noted  that  volume  recombination  and  ionization  would 
increase  the  damping  by  producing  hot  neutrals  and  cold  ions  respectively, 
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Also,  since  most  ion  momentum  is  directed  along  the  z  axis,  elastic 
scattering  of  an  ion  through  anything  other  than  a  small  angle  would 
appear  as  a  damping  effect.   Charge  exchange  would  substitute  a  hot 
neutral  and  a  cold  ion  for  the  original  hot  ion  and  cold  neutral.   Al- 
though the  mean  free  paths  for  charge  exchange  and  ion-neutral  collisions 
are  of  the  same  order  of  magnitude,  the  actual  collision  rate  for  charge 
exchange  would  be  somewhat  less  than  that  for  elastic  collisions,  since 
not  every  ion-neutral  collision  would  result  in  charge  exchange.   But 
charge  exchange  would  certainly  appear  as  another  damping  mechanism. 
All  these  processes  probably  contribute  to  the  inaccuracies  of  the  above 
theory. 

It  should  be  noted  that  the  observed  damping  is  much  less  than  that 
observed  by  Little  and  Jones  in  a  gaseous  discharge  where  the  percentage 
ionization  is  low  and  there  are  many  more  neutrals  present.   5  Also  in 
such  a  discharge  collisions  between  the  ions  and  the  container  walls 
would  cause  additional  damping.   Such  an  effect  is  not  a  factor  in  the 
present  investigation. 

The  primary  effect  of  £  is  to  alter  the  form  of  the  radial  per- 
turbation from"!J"0  to  a  more  rapidly  decreasing  function  of  r.   The  effect 
of  £  on  y  and  z  is  a  complicated  second  order  process;  it  will  not  be 
discussed  further  here. 

Neglecting  the  second-order  effects  (?7>  C"  ,  €.   )  ,  the  equations  4a, 
4b,  and  4c  become  (for  ^A  =0): 

MNVQ  * 


(1*3  -  ?)2  +SL2 


If-  a<o-;» 

MNVr  =  -&-£ 


(i^-  ?)2  +_fl.2 
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MNV„  =  -    ^_ 


iil 


z     (i«0  -9) 

In  this  form  the  equations  clearly  illustrate  the  damping  effect  of  y? . 
If  v  is  now  set  to  zero,  these  equations  can  be  solved  to  give: 


Vz 


.  i    ^>p  N     ojn 


Vr  I     1   >p     o 2 


=  (_±_4HL_) 


Vz      -AP   *r   j!2-02 


V 


z  I     MN«0 

Since  .    6(  Qp    (except  near  cutoff),  it  is  clear  that  most  of  the 
perturbational  momentum  appears  in  the  motion  of  the  ions  in  the  z  direc- 
tion.  If  is  for  this  reason  that  the  waves  are  called  acoustic  waves. 

Theory  for  the  £  ,  JC.    curves  is  derived  from  the  probe  analysis 
mentioned  in  Section  3.   Those  measurements  (which  confirmed  the  pre- 
sence of  only  the  M    =  0  mode)  showed  that  no  density  perturbation  was 
detectable  beyond  15  mm  from  the  plasma  beam  axis.   This  distance  was 
also  found  to  be  independent  of  the  magnetic  field  strength  up  to  2200 
gauss.   This  was  a  remarkably  strong  confirmation  of  the  values  of 
deduced  from  the  data  reduction.   It  also  demonstrates  the  validity  of 
the  basic  dispersion  relation,  since  this  relation  is  far  more  sensitive 
to  small  changes  in  X     than  for  any  other  parameter. 

Fig.  17  is  reproduced  from  the  work  of  Gall  and  Oleson   15  ;  it  can 
be  seen  that  a  value  of  a  =  10  mm  would  be  a  more  reasonable  estimate, 
at  least  from  the  point  of  view  of  electron  temperature  alone.   However 
the  value  of  /C.  -    15  mm  taken  from  the  probe  analysis  gave  a  much  better 
fit  to  the  data  and  was  used  as  the  theoretical  value,  since  it  was 
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independently  confirmed. 

Theory  for  the  6"  ,  q  curves  was  taken  from  the  work  of  Gall  and 
Oleson.  115   Although  their  work  did  not  directly  yield  values  of  O 
a  function  of  the  magnetic  field  strength,  it  did  indicate  that  this 
quantity  is  not  a  sensitive  function  of  the  magnetic  field  strength  in 
the  range  used  in  the  present  investigation  (800  to  2200  gauss).   Fig. 
18  is  reproduced  to  illustrate  this  point.   From  Fig.  18  S     is  seen  to 

be  equal  to  100  m-1  at  B  =  600  gauss.   Assuming  o  is  not  a  function 

p  _   c     _   .303 
of  B0,  £  = 


Bo 
Table  2  summarizes  the  results  of  the  present  investigation  for  the 

low,  medium,  and  high  field  strengths  used. 

Table  2 


quantity  800  gauss       1600  gauss      2200  gauss 

wavelength  (m)  .0458  to  1.19    .0314  to  .341    .0248  to  .145 

i/e  attenuation  dist.      20  to  500       133  to  1110     455  to  2500 
Cm) 

phase  velocity/c  1.25  to  26      1.25  to  9.0     1.2  to  4.4 

group  velocity/c  .09  to  .62      . 1  to  .73       . 1  to  .8 

ion  drift  velocity        90.8  89.4  87.1 

(ms"1) 

effective  ion  collision   342  510  431 

frequency  for  momentum 
loss  (s--1-) 

ion  energy  (ev)  .05  .13  .16 

radial  cutoff  .0152  .0147  .0152 

distance  (m) 

radial  1/e  profile        .0100  .0105  .0105 

parameter  (m-J-) 
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Since  the  values  of  y  and  z  vs.  x  cover  such  wide  ranges,  a  tabula- 
tion of  these  quantities  vs.  x  is  not  included  at  this  point.  Represen- 
tative values   can  be  read   from  the  curves   presented   in   Section  4. 

Table    3 

Parameter  Measured  Value  Deduced  Value 

w|  95.2  m/s 

?  421  s~l 

T  none 

a  15  mm 

S  100  m"1 


85.0 

m/s 

480    ! 

B"1 

.092 

ev 

15.5 

mm 

102  m"1 

99 


Section  6 
CONCLUSIONS 

In  conclusion  it  may  be  stated  that  for  the  propagation  of  ion 
acoustic  waves  where  Landau  damping  is  not  present  the  new  dispersion 
relation  derived  in  Section  2  is  capable  of  predicting  the  propagation 
characteristics  of  these  waves  quite  accurately.   The  greatest  error 
present  in  the  dispersion  and  attenuation  curves'  (that  is  the  greatest 
departure  between  data  and  theory)  is  6%.   It  was,  however,  much  more 
difficult  to  arrive  at  satisfactory  theories  for  the  5  adjustable  para- 
meters.  The  close  agreement  in  the  dispersion  and  attenuation  curves, 
however,  is  believed  to  result  from  the  unambiguous  boundary  condition 
specified  and  from  the  use  of  photomultiplier  detection  (rather  than 
probe  detection  with  its  inherent  errors). 

Where  damping  processes  other  than  ion-neutral  collisions  occur, 
the  effective  ion  collision  rate  for  momentum  loss  may  be  redefined  to 
include  all  such  processes  (such  as  the  collision  rate  of  ions  with  the 
walls  of  the  container).   For  this  case  V?  would  be  the  sum  of  the  in- 
dividual collision  rates  for  all  processes.   Since  negligible  momentum 
is  carried  by  the  electrons,  electron  momentum  loss  mechanisms  may  be 
ignored  without  appreciable  error. 

The  externally  accomplished  excitation  and  dispersion  measurements 
described  have  been  a  valuable  diagnostic  tool  for  the  investigation  of 
the  plasma  characteristics.   It  was  by  this  method  that  the  values  of 
W,  p  ,  0~"  ,  X    and  S  were  deduced  from  the  data.   This  method  could  be 
even  more  valuable  where  internal  measurements  (such  as  probe  diagnos- 
tics) are  impractical. 

It  is  recommended  that  future  investigations  be  made  to  test  the 
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theory  at  very   low  frequencies    (1  kc   to    100  kc) .      Due  to  the   frequency 
limitations   of  the  TAB-7   transmitter  employed,   measurements    in  this 
lower  frequency  region  were   not  possible.      Future    investigations    in  this 
region  would   provide   an  additional   test  of   the  theory   and  additional   con- 
firmation of   the  deduced  values  of  the   5   parameters   considered. 

It    is    also   recommended  that   the   noise  spectrum  of    the  plasma  be 
analyzed.      This   would  provide  additional    information  on  the   naturally 
occurring   interior  phenomena. 

Investigation  of   ion  acoustic  wave   propagation   in  other   rare  gases 
(helium  and   neon  for  example)   would   also   provide   additional   checks   of 
the  theory. 
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